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ABSTRACT 

A type la supernova explosion starts in a white dwarf as a laminar deflagration at 
the center of the star and soon several hydrodynamic instabilities (in particular, the 
Rayleigh-Taylor (R-T) instability) begin to act. In previous work (Ghezzi, de Gou¬ 
veia Dal Pino, & Horvath 2001), we addressed the propagation of an initially laminar 
thermonuclear flame in presence of a magnetic field assumed to be dipolar. We were 
able to show that, within the framework of a fractal model for the flame velocity, the 
front is affected by the field through the quenching of the R-T instability growth in 
the direction perpendicular to the field lines. As a consequence, an asymmetry devel¬ 
ops between the magnetic polar and the equatorial axis that gives a prolate shape to 
the burning front. We have here computed numerically the total integrated asymme¬ 
try as the flame front propagates outward through the expanding shells of decreasing 
density of the magnetized white dwarf progenitor, for several chemical compositions, 
and found that a total asymmetry of about 50% is produced between the polar and 
equatorial directions for progenitors with a surface magnetic field B ~ 5 x 10 7 G, and 
a composition 12 C = 0.2 and 16 0 = 0.8 (in this case, the R-T instability saturates at 
scales ~ 20 times the width of the flame front). This asymmetry is in good agreement 
with the inferred asymmetries from spectropolarimetric observations of very young 
supernova remnants, which have recently revealed intrinsic linear polarization inter¬ 
preted as evidence of an asymmetric explosion in several objects, such as SN1999by, 
SN1996X, and SN1997dt. Larger magnetic field strengths will produce even larger 
asymmetries. We have also found that for lighter progenitors (i.e., progenitors with 
smaller concentrations of 16 O and larger concentrations of 12 C) the total asymmetry 
is larger. 
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1 INTRODUCTION 

According to the current modelling, the explosion of a type 
la supernova starts as a thermonuclear deflagration of a 
Chandrasekhar mass white dwarf of carbon-oxygen (C+O) 
or oxygen-neon-magnesium (O+Ne+Mg) compositions. Due 
to the action of several hydrodynamic instabilities, in partic¬ 
ular, the Rayleigh-Taylor (R-T), the initially laminar prop¬ 
agation does not survive as such and the combustion front 
rapidly develops a complex topology. A cellular stationary 
combustion (followed by a turbulent combustion regime) is 
rapidly achieved by the flame and maintained up to the end 
of the so-called flamelet regime when a transition to det¬ 
onation may occur. As in the case of chemical laboratory 
flames (Gostintsev et al 1988), thermonuclear flames prob- 
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ably develop a self-similar behavior due to the action of the 
hydrodynamic instabilities. Under this hypothesis of a self¬ 
similar deflagration regime, the effective flame speed is well 
described by a fractal scaling law, as previously proposed 
by several authors (Woosley 1990, Timmes & Woosley 1992, 
Niemeyer & Hillebrandt 1995, Niemeyer & Woosley 1997, see 
also Ghezzi, de Gouveia Dal Pino & Horvath 2001), namely 

/ r \ (et— 2) 
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Vfrac — Vlam ( I ? (1) 
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where vi a m is the laminar velocity (as obtained, for exam¬ 
ple, in Timmes & Woosley 1992); and Imin are the 

maximum and minimum Rayleigh-Taylor wavelength per¬ 
turbations, respectively (see Chandrasekhar 1961); and d is 
the fractal dimension, which may assume a range of values, 
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2 ^ d < 3 (Woosley 1990; see also Ghezzi, de Gouveia Dal 
Pino & Horvath 2001). 

In a previous work (Ghezzi, de Gouveia Dal Pino & Hor¬ 
vath 2001, hereafter, Paper I), we addressed the propagation 
of an initially laminar thermonuclear flame in presence of a 
magnetic field assumed to be of dipolar geometry. The main 
result of that work was that within the framework of the 
fractal models for the flame velocity in the wrinkled and tur¬ 
bulent regimes, the front is affected by the field that inhibits 
the growth of the R-T instability in the direction perpendic¬ 
ular to the field lines. As a consequence, an asymmetry is 
established between the magnetic polar and the equatorial 
axis that gives a prolate shape to the burning front (see Fig¬ 
ure 1). In that work, an estimate of the intrinsic asymmetry 
as a function of the field and the core density was presented 
and discussed in the context of type la supernova explo¬ 
sions. Nevertheless, several aspects of this problem needed 
to be clarified as, for instance, a computation of the to¬ 
tal integrated asymmetry, as the burning front propagates 
through the outer, density-decreasing shells of the magne¬ 
tized expanding progenitor star, and also its meaning for the 
observation of SNIa remnants. We here present the results 
of one-dimensional numerical calculations of the integrated 
asymmetry through white dwarf progenitors with different 
initial chemical compositions, taking into account the ex¬ 
pansion of the star (which is assumed to be homologous). 


2 INSTANTANEOUS AND CORRECTED 
ASYMMETRY 

According to the model developed in Paper I, in the pres¬ 
ence of a dipolar magnetic field, the ratio of the propagation 
velocity of the burning front in the polar and equatorial di¬ 
rections of the star is given by 


into account the deformation of the flame) may be defined 
by 

A{p) = Vpol{p ~ Ap) . (3) 

Veq(p) 

Since the fractal propagation velocity is a continuous 
and differentiable function of the density, we can expand 
Vp 0 i in a Taylor series around p and keep the zeroth and 
first order terms provided that A p <C p, to yield 
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which may be rearranged to display the instantaneous asym¬ 
metry and a correction term as 


A(p) = A{p) 


1 dvpoi (p) 
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where 
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is the correction term for the asymmetry (see below). 

Given that v po i > v eq and both velocities increase with 
decreasing density *, the difference between v po i and v eq is 
necessarily larger than the difference that they would have 
by imposing the same density ahead for both flames. In other 
words, the corrected asymmetry A(p) must be larger than 
the instantaneous asymmetry A(p). Moreover, the value of 
A c (p) depends on the value of Ap, which in turn, depends 
on the value of the corrected asymmetry A(p) defined in 
Eq.(5). We may then find a solution to the non-linear Eq. 
(5) by discretizing and solving it iteratively 2 
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where v po i and v eq are the velocities at the polar and equa¬ 
tor directions, respectively. This equation shows that the 
asymmetry of the velocity field is controlled by the quotient 
of the magnetic pressure (B 2 /8iv) and the ram pressure of 
the gas (oc pvf arn ). We shall label A(p) as the instantaneous 
asymmetry as a function of the density p for an initially 
spherical flame. In paper I, we have computed this quantity 
for different magnetic field intensities. 

It is clear, however, that A(p) will be increasingly dif¬ 
ferent from the quotient above as time elapses and the flame 
propagates outward through the star, since the flame in both 
directions will be encountering decreasing densities ahead. 
Therefore, a full integration of the asymmetry is needed to 
account for this fact. 

Let us first assume that the star does not expand to 
derive a zero-th order correction. Since the polar flame front 
is faster than the equatorial front, the former will encounter 
lower densities first. If label the difference of densities in 
each direction as Ap = ( p eq — p po i) > 0, the polar front 
faces a density p po i = p — Ap when the equatorial front is at 
pe q = p. To this order, a corrected asymmetry (which takes 


An{p) = A„_ 1 (p) - ( Ap„_!, (7) 

V v eq ) n -1 

with A 0 = A(p), and A p a = 0. Whenever Ap ~ p, higher or¬ 
der terms become important, although p Ap is automat¬ 
ically satisfied. Eq.(7) converges provided that the Cauchy 
condition 

\A n (p) - A n _i(p)| << e (8) 

is satisfied for an arbitrarily small e. We shall see below 
how to integrate the corrected asymmetry through the star 
so as to estimate the accumulated effect along the whole 
propagation of the non-spherical flame. 


2.1 Correction Term of the Asymmetry 

The correction term for the asymmetry (Eq. 6), contains a 
derivative of the polar speed with respect to the density, and 
therefore, we need to know the dependence of the laminar 

1 The fractal velocity, Vf rac , increases with the density, in Eq. 
(1), because the minimum scale l m in decreases faster with density 
than the laminar velocity, vi am (see Paper I). 

2 It can be checked that the Lipschitz condition is satisfied by 
the function A(Ap). 
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speed with the density and composition of the fuel in order 
to calculate the correction term. In this subsection we cal¬ 
culate the correction term for C+O progenitors. Using the 
interpolation formulae given by Woosley (1986); Timmes & 
Woosley (1992); and Arnett (1996), we find 


Vlam = 92 x 10 5 


/ \ 0.805 

( P \ 

'X( 12 C)‘ 

\2 x 10 9 / 
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cms \ 

which is approximately 3 valid in the density range 
0.01 ^ pg ^ 10. In the particular case in which d = 2.5 4 , 
the effective polar speed v po i (see Eq. 1) does not depend on 
the laminar speed viam, so that we could say that v po i is to 
some extent “independent of the microphysics” 5 
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And, in general, the correction term is given by 
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where, the value of A p must be calculated numerically (see 
the next section). 


v P oi(pA=0.282p 0 u s ( 9Lm 2 a ; 5p )°' 5 cms" 1 . (10) 
V Pu J 

Here, dp = p u — pb is the density jump at the flame, i.e., the 
difference between burned ( pb ) and unburned ( p u ) material 
densities. We note that this is, in general, different from 
the the density difference between the polar and equatorial 
fronts A p = p eq — p po i (see Eqs. 3 - 7). Thus 
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Using Eqs. (1) and (9), it is possible to obtain the po¬ 
lar fractal speed, v po i, for any C+O composition with an 
arbitrary fractal dimension d 


3 INTEGRATED ASYMMETRY 

In order to compute the total asymmetry for a given set 
of initial conditions for the progenitor, we must integrate 
the following equations along both the equatorial and polar 
directions, respectively 


dr, 




(14) 


dr po i 
dt 


= V po l(p,t) 


(15) 


where v po i(p) = A(p) v eq {p), and v eq is given by: 


Veq — Vl arn 


le 


(d-2) 


(16) 


where l eq = l m i n at the magnetic equator of the star (Paper 

I): 


Vpoi = 14.1683 (0.208752) d X( 12 C)°' 889 x (12) 


3 We should note that equation (9) is only approximate and it is 
used for simplicity reasons, more realistic results could probably 
be obtained directly interpolating the velocities given in table 
(3) of Timmes & Woosley 1992, and will be adressed in future 
calculations. 

4 This is actually the maximum fractal dimension estimated in 
Paper I. 

5 In fact, recent numerical simulations of supernovae have shown 
that the hydrodynamics of the explosion is independent of the mi¬ 

crophysics of the flame (Gamezo et al. 2002). In paper I, we found 
that the minimum relevant hydrodynamic scale for the flame was 

lmin ~ 10 5 cm (for p rsj 10 9 gem -3 ), and the flame evolution was 

’’independent of the microphysics” since lmin Sf almost all 
the way, where 8f is the flame width. We here stop the numerical 

integration when lmin ~ 8f or lmin ~ 20 x 5f (see below). In the 

particular case that d = 2.5, using the fractal scaling given by eq. 
(1) we recover the independence with microphysics which was ob¬ 
tained through dimensional arguments by Khokhlov (1995), and 
is also consistent with laboratory results of Taylor (see Kull 1991). 
In this case, the flame does not depend directly on quantities, such 
as, the laminar velocity or the thermometric conductivity. 


(17) 

The equations above can be integrated simultaneously 
in one dimension outward through the star along with Eqs. 
(7) and (8) to give the integrated asymmetry A to t(p) = 
r po i(p)/r eq (p), up to a shell of density p. The integration 
of Eq. (14) is straightforward, while Eq. (15) is non-linear 
since the value of A depends on the actual position of the 
polar front r po i. It is solved using the same iterative algo¬ 
rithm that gives the corrected asymmetry described in the 
previous section. 

Until now, we have neglected the fact that the flame 
is actually propagating in an expanding star and have thus 
partially decoupled the equations of the burning front from 
the hydrodynamical motions of the full star. Since actual 
deflagrations are subsonic, the star will pre-expand ahead of 
the flame, and this will, in turn, affect the burning itself (see 
Figure 1). 

To a very good approximation, we can assume an ho¬ 
mologous expansion of the white dwarf and this allows a 
simple description of the expansion as a function of time, in 
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which the radius of the n-th zone at a given time t evolves 
according to 

r„(t) = a(t) r° , (18) 

where r° gives the distance to the center in the initial condi¬ 
tion (i.e., in hydrostatic equilibrium), and a(t) is the homol¬ 
ogy factor satisfying a(0) = 1. The density will then evolve 
according to 

p n (t) = p„/a 3 {t), (19) 

where p ° is the density at the n-th zone in hydrostatic equi¬ 
librium, so that in the comoving frame of each zone the 
density will decrease with time. 6 

It is easy to see that the stellar expansion should pro¬ 
duce a slower flame front than in the stationary case. This 
is due to the fact that when the deflagration reaches the 
radius r n , the density there will be lower than it is with¬ 
out expansion. Therefore, the propagation velocity will be 
modified accordingly. The distance between zones is now 
dr' eq = a(t) dr eq , and dr po i = a(t) dr po i , for each of the prop¬ 
agation directions. The n-th zones evolve according to 

r '?n = Ve( ln — 1 + d-r e q n _ 1 , 

and 

r poi n = T po i n _ 1 + dr po i n _ 1 , 

respectively. 

We have computed the total asymmetry by taking into 
account the stellar expansion effect upon the density, for dif¬ 
ferent initial chemical compositions of the white dwarf pro¬ 
genitor. We have integrated Eqs. (14) and (15) together with 
Eqs. (18) and (19), using a subroutine to interpolate step by 
step the physical properties of the deflagration front, i.e., the 
laminar velocity and the density jump, and the value of the 
magnetic field. We have employed an analytical function for 
a(t), namely a(t) = 0.25 1 + 1 obtained from a numerical fit 
to the homologous expansion solutions of Goldreich & We¬ 
ber (1980). The magnetic field strength has been assumed 
to decrease approximately linearly with the radial coordi¬ 
nate (see Wendell et ah 1987, and Paper I), and the initial 
density and gravity of the progenitor are depicted in Figures 
2 and 3, respectively, as functions of the stellar radius. To 
compute the variation of the magnetic field in the comoving 
frame of the flame, we have assumed magnetic flux conserva¬ 
tion. The computed progenitor models were resolved in 10 3 
zones, giving a physical resolution ~ 1.5 km for a star with 
a radius R= 10 5 km. The integration was stopped when the 
polar combustion front reached a shell with density p ~ 10 7 
g cm' 3 for which the transition to detonation is believed to 
occur (Khokhlov 1995). This is a first step towards a more 

6 We note that some deviation of a perfect homologous expan¬ 
sion should be expected mainly due to the fact that, as the flame 
is accelerated downstream under the action of hydrodynamical 
instabilities, shock and pressure waves may travel ahead of the 
flame and disturb locally the stellar structure. Nonetheless, since 
the deflagration evolves under approximate isobaric conditions, 
the actual expansion is expected to be nearly close to the homol¬ 
ogy assumed here. 



Figure 1. Schematic representation of the propagation of the 
front inside a white dwarf. Here B indicates the direction of the 
dipolar magnetic field, V e xp is the star’s expansion velocity, V po ; 
and ~V e q are the fractal polar and equatorial velocities, and r po i 
and r eq are the polar and equatorial radii, respectively, as de¬ 
scribed in the text. We note that the combustion is isobaric to a 
very good approximation, and the density jump in the flame is 
~ 10 —1 times the density value. 



Figure 2. Initial density of the white dwarf progenitor as a func¬ 
tion of the stellar radius. 

complete calculation, but nevertheless contains all the es¬ 
sential ingredients necessary to address the total asymmetry 
A’tot, namely 

A 'tot = r' po i/r' e q , (20) 

where, generally speaking, r' eq > r eq and r pol > r po i. 


4 RESULTS 

The results of the calculations are shown in Figures 4 to 8. 
The instantaneous asymmetry near the stellar surface (for 
a density p — 5 x 10' gem -3 ) is shown in Figure 4 for sev¬ 
eral progenitors with C-0 compositions as a function of the 
surface magnetic field. 

Integrated asymmetries are depicted in Figure 5 for 
three different compositions (see captions) and a maximum 
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Figure 3. Initial absolute value of the gravity of the white dwarf 
progenitor as a function of the stellar radius calculated for a pro¬ 
genitor with an electron fraction Y e = 0.5. 
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Figure 4. Instantaneous asymmetry for progenitors with dif¬ 
ferent initial compositions [ (a) A ( 12 C) = 0.2, X( 117 O) = 0.8; 
(b) X{ 12 C) = 0.5, X( le O) = 0.5; and (c) X( 12 C) = 1.0, 
X( le O) = 0.0] as a function of the magnetic field at a density 
p = 5 X 10 7 gem -3 . 


surface magnetic field B= 10 9 G. It is apparent from Figs. 4 
and 5 that heavier progenitors develop larger instantaneous 
asymmetries, but the lighter ones achieve the largest inte¬ 
grated asymmetries. Hence, it is important to understand 
the reasons behind such behavior for future applications. 

The key aspect for this behavior may be traced back to 
saturation effects of the flame as follows. The value of l m in 
(see Eq. 17) indicates that it must decrease as the density 
decrease, but it clearly cannot be arbitrarily small. In fact, 
there is a minimum value associated with the width of the 
flame 5/, which we will call the saturation scale, l sa t- 

The asymmetry between the polar and equatorial direc¬ 
tions begins to decrease when the instabilities in the polar 
front reach the saturation limit l 3a t- The polar front will then 


B s = 10 9 G 



Figure 5. Integrated asymmetries for the progenitors (a) 
X( 12 C) = 0.2, X( le O) = 0.8; (b) X( 12 C) = 0.5, X( 16 0) = 0.5; 
(c) X( 12 C) = 1.0, X( 16 0) = 0.0. A surface field B= 10 9 G is 
assumed for all cases (which corresponds to a core B~ 10 10 G), 
and a saturation scale of 20 5/. Taking lower values for the satura¬ 
tion scale higher integrated asymmetries are obtained. A general 
trend of the asymmetry effect is seen here, i.e., lighter progeni¬ 
tors show higher integrated asymmetry (while heavier progenitors 
have higher instantaneous asymmetry, for a given density and 
magnetic field, as shown in Fig. 4). This graph can be also inter¬ 
preted as a display of the evolution history of the asymmetry. In 
the beginning, when the flame is still propagating in a high den¬ 
sity medium, the asymmetry is very low. Later on, near the end 
of the deflagration regime, the asymmetry largely increases until 
the front encounters very low densities and the flame vanishes. 


stop accelerating, while the equatorial front will be still ac¬ 
celerating. This will give place to a symmetrization phase 
that will last until the flame is quenched, or a transition 
to detonation occurs. We cannot predict exactly this sym¬ 
metrization effect here, but if a transition to detonation oc¬ 
curs at densities ~ 10' g cm -3 , then the errors in our present 
calculation must be small. On the other hand, if no transi¬ 
tion to detonation takes place, then the uncertainties in our 
estimates will depend upon the density at which the flame 
quenches and will be as large as small the density is. Due to 
the same reason, we predict smaller symmetrization effects 
on lighter progenitors, i.e., progenitors with V( 12 C) = 0.5 
and V( 16 0) = 0.5 will reach the saturation limit later on at 
lower densities, and will, therefore, develop larger integrated 
asymmetries (see Fig. 5). 

It is known that if l sa t = lmin ~ <5/ thermodiffusive ef¬ 
fects must be taken into account and, in fact, in most of the 
laboratory flames lsat -5 20 5/ is observed (see Bychkov et 
al. 1999). Since no definitive answer to this question exists, 
we have calculated the integrated asymmetry for different 
values of the saturation scale between these two extremes 
(i.e., from l sat = 5f to 20 5/) in Figure 6, for the progenitor 
(a) of Figure 5. In Figure 7, we have fixed the saturation 
scale at l sa t = 20 5/, in order to compare the integrated 
and the instanteneous asymmetry variations with the sur¬ 
face magnetic field for the same progenitor. Analysing these 
figures, we find that the integrated asymmetry is enormous 
if lsat ~ Sf, and can be larger than 1000% if the magnetic 
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B s = 5 x 10 7 Gauss, 12 C = 0.2, “O = 0.8 



Figure 6. Integrated asymmetry vs saturation scale l sa t in units 
of the flame width 5/, for the progenitor (a) of Figure 5 with 
chemical composition X( 12 C) = 0.2, X( 16 0) = 0.8. A surface 
field B= 5 X 10' G is assumed. 


12 C = 0.2, 16 O = 0.8 



Figure 7. Circles - Integrated asymmetry vs surface mag¬ 
netic field strength for the progenitor with chemical composition 
X( 12 C) = 0.2, X( 16 0) = 0.8, and a saturation scale l sa t = 20 5/; 
Full line - seven-degree polinomial curve that fits the numerical 
results; Dashed line - instantaneous asymmetry for the same pro¬ 
genitor at a (surface) density p ~ 10 7 g crn " . 

field is > 10 9 G. When l sa t ~ 20 5/ is imposed, the asym¬ 
metries are smaller but can be still very large, for instance, 
the progenitor shows a maximum asymmetry of the order 
of 350% for B = 10 9 G (Fig. 7). When the field is reduced 
to 5 x 10 7 G, the asymmetry drops. In this case, it presents 
a total asymmetry of ~ 50%, for l sa t ~ 20 5/ (Figure 8), 
which is in the ballpark of what polarimetric observations 
suggest (~ 20 % or less, see the next section; Leonard et al. 
2000, Wang et al. 1997, Howell et al. 2001). 


B s = 5 x 10 7 Gauss, 12 C = 0.2, 16 Q = 0.8 



Figure 8. Integrated asymmetries for progenitor (a) AT 12 C) = 
0.2, X( 16 0) = 0.8, with a surface field B = 5 X 10 7 G, and sat¬ 
uration scales l a at = 5/ (dashed line), and 20 5/ (full line). This 
represents the set of parameters that better fits observational data 
from SN la spectropolarimetric studies and experimental knowl¬ 
edge from laboratory flames. We have stopped the calculations 
when the flame vanishes. At this time, expansion can start sym¬ 
metrizing the flame, although more efficient symmetrization must 
be expected during the coasting phase of the remnant due to en¬ 
hanced turbulent diffusion. 

seems to be evidence of asymmetric explosion (see Branch 
et al. 2001; Leonard et al. 2000; Howell et al. 2001). In the 
present work, we have investigated the effects of the mag¬ 
netic field of the white dwarf (WD) progenitor on the prop¬ 
agation of the burning front of thermonuclear supernovae, 
and found that an asymmetry develops that could, in princi¬ 
ple, help to explain observed asymmetries in these systems. 

The magnetic field strengths inferred for isolated WDs 
are between 3 x 10 4 G and 10 9 G, while for WDs in AM Her 
binaries, they are in the range 10' G to 2 x 10 8 G (Wick- 
masinghe & Ferrario 2000) and could be further amplified 
by compression during the accretion phase (Cumming 2002). 
Magnetic fields had not been taken into account in previous 
studies of thermonuclear supernovae because their pressure 
is found to be much smaller than the gas pressure. Nonethe¬ 
less, our studies have revealed that the magnetic fields act 
upon the hydrodynamical instabilities that develop in the 
flame front (where the magnetic pressure is larger than the 
ram pressure; see Eq. 2) and quench their growth in the di¬ 
rection perpendicular to the field lines. As a consequence, an 
asymmetry develops between the magnetic polar and equa¬ 
torial axis that gives a prolate shape to the burning front. 

We should emphasize that our analysis of the Rayleigh- 
Taylor mode interactions at the flame front has been based 
on the model described by Zeldovich et al. (1966, 1980) 7 , 
and its results are in qualitative agreement with 2-D and 3- 


5 DISCUSSION AND CONCLUSIONS 

It has been generally assumed that supernovae explosions 
are spherical. However, recent polarimetric studies of several 
supernovae have revealed intrinsic linear polarization that 


7 The Zeldovich et al.’s model for the mode-mode interactions 
of burned cells explains the saturation of the exponential growth 
of the hydrodynamic instabilities and gives the lower cut-off as¬ 
sumed for the instability which depends mainly on the saturation 
scale as explained above (see Eq. 17 above and also Eq. 4 of Paper 
!)• 
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D numerical studies of the Rayleigh-Taylor (R-T) instabil¬ 
ity development in fully turbulent magnetized layers without 
combustion (see, Jun, Norman & Stone 1995, and references 
therein), and with numerical simulations of thermonuclear 
flames in 2D and 3D (see Khokhlov 1995). Jun et al.’s re¬ 
sults, actually reveal a tendency, in the nonlinear regime, for 
an increase of the growth rate of the Rayleigh-Taylor insta¬ 
bilities in the direction parallel to the magnetic field lines 
(the polar direction), which would amplify the asymmetry 
effect here investigated, and an inhibition of the growth rate 
in the direction perpendicular to the magnetic field lines 
(equatorial direction) like in the present investigation. 

The calculations performed in the former sections show 
that the integrated asymmetry in the velocity fields may be 
very large if the surface magnetic fields are near a maxi¬ 
mum value (~ 10 9 G). Observed asymmetries are much less 
extreme. At least three alternatives to solve this mismatch 
come to mind. The first is that white dwarfs in binary sys¬ 
tems, which are thought to be the progenitors of SNIa super¬ 
novae, never have magnetic fields of such an intensity. Since 
pre-supernova systems samples are highly incomplete, it is 
difficult to address the likelihood of this conjecture (see be¬ 
low). The second possibility is that asymmetries do not de¬ 
velop that much because the saturation scale l sa t ^>20 5/ in 
these systems. This possibility cannot be totally disregarded 
since the saturation scale may be substantially different in 
these systems from that obtained in laboratory deflagration 
fronts. The third way out would be that even if the explosion 
itself could be largely asymmetric, the observed remnant 
could be efficiently symmetrized in the very early phases 
of the expansion. Nevertheless, it would be worth to search 
for a subset of very asymmetric remnants, perhaps not even 
identified as such because a jet-like morphology has never 
been expected. 

On the other hand, if the surface magnetic field is re¬ 
duced to more moderate values (< 5 x 10' G), we have 
found that the total asymmetry drops (for saturation scales 
~ 20 5/) to a value of ~ 50% which is much more in agree¬ 
ment with the asymmetries inferred from present observa¬ 
tions. For instance, for the type la supernova 1999by, the 
first clear example of an object of that class with asymmet¬ 
ric structure, an intrinsic polarization of ~ 0.8% was ob¬ 
served that was interpreted as produced by radiation from 
an oblate or prolate distribution of scattering electrons with 
an asymmetry between 17% and 50%, depending on the in¬ 
clination of the object relative to the line of sight (Howell et 
al. 2001). Another important aspect to be noted is that in 
SN1999by, Si lines stronger than in normal type la SNs have 
been detected. This could be explained, in the framework of 
our model, as a result of an incomplete burning in the equa¬ 
tor direction followed by a smaller production of Ni which, 
in turn, could explain the observed redder light curve. Yet 
another possibility may be that the prolate flame produces 
a transition to detonation first on the polar direction, thus 
disrupting the star before the equatorial direction had the 
chance to reach the requirements for detonation transition, 
therefore resulting in an incomplete burning and smaller Ni 
production. 

So far, only ~ 15% of the type la supernovae have been 
examined through polarimetric observations and they do not 
present, in general, high polarization degrees. The model 
here described, could, in principle, explain the asymmetry 


inferred from these spectropolarimetric observations, if the 
fraction of strongly magnetized white dwarfs (WD) is not 
too large. In fact, it seems that only 5% of the isolated white 
dwarfs are magnetized, but this figure increases to 25% in 
binary systems. 

Very recently, Gumming (2002) has suggested that due 
to an effect of screening and magnetic held amplification in 
rapidly accreting systems, an asymmetry (due to magnetic 
held effects) can be more likely observed in super-soft X-ray 
binaries. This constitutes a new (and separate) observational 
constraint to the study and detection of the asymmetry ef¬ 
fect presented here. In fact, the search of asymmetric SNIas 
in super-soft X-ray binaries could be an easier task than to 
survey a full sample of very young type la supernova rem¬ 
nants through spectropolarimetric observations. 

In the model proposed here, if higher expansion coeffi¬ 
cients, a(t), were considered, then either a higher magnetic 
held strength or a lower saturation scale would be needed 
in order to produce the same value for the total asymmetry. 
In other words, we are not able to establish completely a set 
of parameters for the model, and a statistically significant 
sample of asymmetric type la SNs would be required for that 
purpose. It must be also pointed out that in a real explosion, 
a(t) should assume a distinct functional dependence on the 
pole and the equator, and this would enhance the asymme¬ 
try effect. In the model presented here, there is, in fact, a 
spherical symmetric remnant surrounding a prolate flame. 

We have assumed in the model a constant fractal dimen¬ 
sion all over the explosion. This approximation is justihed by 
the very fast growth of the instabilities in the thermonuclear 
front, so that the turbulent regime is achieved well before 
an appreciable expansion of the star. 

Finally, we should mention that other models have been 
suggested in the literature to explain the asymmetries in 
type la supernovae explosions, although all of them have 
their problems (see e.g., Howell et al 2001). The model here 
investigated also presents limitations that could, at least in 
part, be solved with the help of future multidimensional nu¬ 
merical calculations and larger samples of observed asym¬ 
metric astrophysical objects. Laboratory experiments of the 
deflagration of mixtures of gases in the presence of magnetic 
fields could also provide alternative tests of the asymmetry 
effect in the flame front and also of the fractal model exam¬ 
ined. The configuration of these experiments must be settled 
carefully in order to ensure the growth of the R-T instability 
over a wide spectrum of perturbations, and the geometry of 
the apparatus must be chosen in such a way that its min¬ 
imum size is larger than the minimum wavelength of the 
perturbations (Ghezzi 2002). 8 
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8 We are of course supposing that the fractal model can explain 
the turbulent regime of a supernova flame. Previous laboratory 
experiments (Gostintesv et al 1988) and numerical simulations of 
fractal flames in the wrinkled regime (Filyand et al. 1994; Blin¬ 
nikov et al. 1996) seem to support this hypothesis. 



8 C. R. Ghezzi 
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